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We consider a bilayer version of the extended t-J model, with a view to computing 
the form of certain experimentally observable properties. Using the slave-boson decom- 
position, we show at the mean-field level that in the bilayer system the existence of 
in-plane d-wave singlet pairing excludes any interplane singlet order for reasonable val- 
ues of the interplane superexchange parameter. Restricting the analysis to the regime of 
no interplane singlet pairing, we deduce parameter sets reproducing the Fermi surfaces 
of YBCO- and BSCCO-like bilayer systems. From these we calculate the form of the 
dynamic susceptibility x(ct) ^) in both systems, and of the anomalies in frequency and 
linewidth of selected phonon modes in YBCO. We compare the results with experiment, 
and discuss the features which differ from the single-layer case. 
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1 Introduction 



Bilayer models for high-Tc superconductors have become popular recently, not only due 
to the obvious feature that the best-characterised material has a bilayer structural unit, 
but also since the suggestion by Millis and Monien [||] that the single- and bilayer systems 
may have fundamentally different magnetic excitations. These differences may include 
the existence or absence of the spin gap p, 0], which has become a crucial experimental 
feature and theoretical benchmark. 

In this work we will pursue an approach based on the t-J model |^ analysed by the 
slave-boson decomposition j^, 0, 0. In the mean- field approximation, this approach 
has been studied intensively fl^, |ll], |l3], |l5], |l3 with a view to describing the 



magnetic properties of the high-T^ materials. As an example of such investigations, the 
extended t-J model of a single Cu02 layer, has been shown |Tl|, [l^ to give a good account 
at the RPA level of many features of the spin excitations in the monolayer La2_xSra.Cu04 
(LSCO) system, and some of those in the bilayer YBa2Cu307_5 (YBCO) case. These 
authors considered the (in)commensuration of magnetic fluctuations, the frequency de- 
pendence of the dynamic susceptibility and the different temperature dependences of 
the shift and rate of nuclear magnetic resonance between high- and low-doping regions. 
Such models also contain a consistent description of the "spin-gap" phenomenon , first 
noted by Yasuoka |0. Taking fluctuations about the mean-field solution into account by 
a gauge- field approach [0 , it is possible to gain an understanding of various transport 
properties, including the temperature dependence of the resistivity, thermopower and 
Hall coefficient. Again within the same framework, it has recently been shown [19| that 
one may obtain a good account of several lattice-related features of the electronic sys- 
tem, such as phonon anomalies and an isotope effect, by considering the coupling which 
arises naturally between phonon modes of the layer and the spin sector. 

In view of this degree of correspondence to experiment, we consider the extended 
t-J model to be one of the leading candidates for a framework in which to construct 
a coherent understanding of the many and complex features of the high-Tc problem, 
and thus that the detailed computation of physical quantities within it is a valuable 
exercise. Here we wish to apply the ideas of the single-layer studies and [|l^ to 



a consistent bilayer model, and thus to elucidate the successes and limitations of the 
mean-field approach to spin-dependent, microscopic properties. We draw attention to 
a brief study [0 which estabhshed that the bilayer dynamic susceptibility has the 
experimentally observed periodicity. Ubbens and Lee [l^ have also considered static 
magnetic properties in a bilayer, nearest-neighbour t-J model, including gauge-field 
effects, and argue for an anisotropic (or "extended") s-wave pairing state of coupled 
intra-and interplane singlet order, as well as for enhanced spin-gap formation. Similar 



results for the spin gap were obtained in Ref. pO |. 

The question of the gap symmetry has remained a topic of much debate: while it 
is known in the single-layer model that the low-energy state is d-symmetric |^, |^, a 
variety of reports find somewhat different behaviour in the presence of coupling to a 



second layer. In addition to the result of Ref. [16|, the weak-coupling antiferromagnetic 
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spin fluctuation theory, which favours d-symmetry in a single plane, has also been shown 



to favour extended s-symmetry in a bilayer pT|. A recent numerical analysis by Eder 



et al. P2| gives the evolution of the in-plane cZ-wave state as interplane coupling is 



increased, and demonstrates in a small cluster that there is an abrupt crossover to a 
predominantly interplane ordered state at some value of the coupling. We will discuss 
these results in comparison with our own in section 2. 

The bilayer spin system has also been studied extensively by Millis and coworkers 
^ with the aim of elucidating the quantum nature of the ground state and 
explaining the contrasting NMR results in single- and bilayer materials. From spin-wave 
and scaling theories it has been shown that interlayer coupling acts to provide a 
strong enhancement of the spin gap state (by suppression of competing magnetic order), 
which provides an explanation of the differences without the requirement that the single- 
and double-layer systems have completely different ground states. These considerations 
have been cast in a more microscopic form to describe electrons on a bilayer with a 
weak interplane antiferromagnetic interaction ||20[: the system may vary between a 
Fermi liquid state, which is proposed to describe the doping regime of materials with 
no spin gap, and a spin liquid picture of fermions interacting via a gauge field, which 
shows spin-gap features. 

The structure of this paper is as follows. In section 2 we consider the mean-field 
equations of the bilayer problem, deduce the nature of the gap and solve for BSCCO- and 
YBCO-like systems. In section 3 we calculate the dynamic susceptibility for both types 
of system, and illustrate the dependence of the results on the wavevector component g^. 
In section 4 we compute the superconductive anomalies in c-axis phonon modes of in- 
plane oxygen atoms for the YBCO system, showing the contrast between modes of even 
{g) and odd {u) symmetry. Section 5 contains a summary and concluding discussion. 



2 Mean-Field Solution 



We consider a three-dimensional system consisting of a periodic but uncoupled stack of 
coupled Cu02 bilayers, in which the intra- and interplane hopping (t) and superexchange 
(J) interactions are as shown in Fig. 1. Denoting the c-axis dimension of the unit cell 
by c and the separation of the planes of the bilayer by the ratio Tc = d/c will appear 
in a phase factor arising from the bilayer spacing. The Hamiltonian for a system of Nz 
bilayer units is 



H 




n=l 



ija (ij) 



It 2 I 2t 1 



(1) 



where / is the layer index in each unit cell. The first line contains the t-J model of a 
single Cu02 layer, where tij denotes the extended transfer integrals, Jij the superex- 
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change interaction, which is assumed to be finite only between nearest neighbours, and 
the Hilbert space excludes double occupancy of the quasiparticles ai„. The second line 
contains the coupling between the layers, which may be by both hopping and superex- 
change. Within the plane, the nearest-neighbour hopping term is taken to be t = 4J, 
while second- and third-neighbour terms will be chosen below in order to match the 



Fermi surface shape of the physical systems |Tl|]. We will consider initially a range of 
values for the interlayer coupling parameters and Jj_. 

The assumption of coherent, single-particle hopping processes between the planes 



of the bilayer is counter to the interlayer pair tunnelling theories |25|, which demands 
that this vanish by orthogonality. A clear experimental indication of whether interlayer 
transfer is coherent or incoherent should be given by photoemission experiments, which 
will see either two Fermi surfaces or one. At the moment there is some controversy, and 



inconsistency, surrounding the interpretation of results obtained for both BSCCO |26 



P7[ and YBCO [^, ^ systems, part of which centres on the issue of the number of true 
band-crossings observed where the bonding and antibonding bands of the coherently- 
coupled bilayer are expected to be far apart in the Brillouin zone. Here we avoid further 
discussion of this topic and proceed to elucidate the properties of the bilayer system with 
coherent hopping, but will refer below to the 2-band interpretations of the experimental 
data. 



In the slave-boson decomposition, the operator ajs is represented as ajs = fjsb], 
where fjs is a fermion (spinon) carrying the spin degrees of freedom Sj, and bj a bosonic 
holon carrying the charge. In the mean-field approximation, the spin and charge behave 
independently, and all 4-operator terms may be written expressions quadratic in the 
operators with the mean-field order parameters (finite expectation values of possible 
2-operator combinations) as coefficient, while terms quadratic in these order parameters 
are also generated. Here we take as possible in-plane order parameters 

i) Xi,i+r = EaifilfUra)- uulform RVB of spinous 

ii) A,,+, = j^.in^fUn - fliflr^Y singlet RVB 

iii) -Bi,j+r = {b\J)%ra)'- uulform RVB of holons 
and as interplane order parameters the combinations 

i) X»j = Eaiflafja) = J2a{fiafja)- mtcrplane uniform RVB of spinous 

ii) = W'^fh - fhfp = T^if'^fh - fnfhy interplane singlet 

iii) B'^j = (bj^b^l^) = {b]^b^l): interplane uniform RVB of holons. 

For historical reasons the terminology "uniform RVB" is taken to signify the hop- 
ping order parameter for spin or charge, whose finite value signifies diagonal coherence 
of the motion. As will be shown below, for the spin order parameters it is sufficient to 
consider only those between nearest-neighbour sites, so that r = ±x or ±y for x and A, 
and r = ±rcZ for x' a^nd A'. In the single-layer problem, the free energy is minimised 
when Xi,i+T is uniform, but when Aj j-t^ = — Aj j-ty [^, §|, so that in reciprocal space the 
gap function A^ has cZ-wave symmetry. 
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With these definitions, the mean-field approximation to the Hamiltonian (Q) in re- 
ciprocal space can be represented as 



k k 



(2) 



where Hq contains terms quadratic in the order parameters, and b) 

ffc = {fkpf-kiyfk\jf-ki) the bosou and fermion state vectors, 
fermion matrices are given by 



bl\bl^) and 
The boson and 



M 



bk 



, , IfO 



(3) 



and 



A. 



M 



Afe 




Dl 


Ml 




Dk 


-a 


Afc 




Mk 


Afc 


6 / 



(4) 



with the following definitions. The band energies for the quasiparticles are = ~\tfk + 
^'fk + ^k + jJXk] — /i and Uk = —[tbk + t'bk + ^bk\ ~ -^j which the in-plane hopping 



terms are t 



fk 



J2k' ^' 



k'-k\^k 



,6fc/) and f^, = Ek'^^t' 



M'Jk'a). with tl = ErJ'-e^-^ 



and /i and A are fermion and boson chemical potentials. The (i-symmetric gap is given 
by Afc = ^^JA(coskx — cos ky). The effective uniform interplane order parameter is 
Mk = t±_kB'k + |J±Xfc! where the k dependence of t^k will be discussed in detail below, 
and the interplane gap parameter is = ^^J±A',. 

The 4x4 unitary matrix which diagonalises the Hamiltonian of the fermionic part 
may be written as 



Tfc"^ 



/ cos 6^ sin O'l 
1 — sin 9'^ cosh 6*^ 

71 e^-^cos^^ e^'^sin^+ 

V -e^-^sin^^ e*<^cos^^ 



— e 



cos6'^ 
— sin^i 



'cos^^ -e"^"^sine^ \ 



' cos 01 



sin 9, 
cosh6't 



(5) 



where 



cos 29^ 



a ± \Mk 



and 



sin 29^ 



Afc ± \Dk 
Et 



(6) 



The phase cj) = —qzTc is chosen to cancel the gz-dependent phase of the interplane order 
parameters, and will be the same for both and the modulus in IQfcl thus denotes 
the absence of this phase, but does not prevent the order parameter from being negative. 
It is the transformation (^) which will be used extensively in the sections to follow, when 
the boson operators ffc are replaced in the physical quantities to be computed by the 
quasiparticle operators 7fc = Tfcffc. 
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The fermion part of the Hamiltonian then takes the diagonal form 





M 



( -E^ 




V 



\ 



-E, 

E, J 



(7) 



in which the eigenvalues are given by the combinations 



Et 



(a±|Mfc|)' + (Afe±|D, 



1/2 



(8) 



E'f^ corresponds to the bonding band of the bilayer system, and i?^ to the antibonding 
band. 

The 2x2 bosonic matrix is trivially diagonalised by the transformation 



to new quasihole operators f3k, in terms of which the diagonal boson matrix is 



M 



bk 



-UJi- 



where the eigenenergies are given by 



(9) 
(10) 



(11) 



(12) 



The system is required to satisfy eight mean-field equations, three for the in-plane 
order parameters, three for the interplane parameters and two for the carrier number 
5. These equations may be obtained by minimising the free energy with respect to 
each of the parameters in turn, or by substituting the quasiparticle operators into the 
reciprocal-space definitions of the order parameters and transforming back to real space. 
One obtains for the in-plane parameters 



-ik.r 



2N 51 ' 



V2N 



COS 26*^ tanh ^pE^ + cos 26^^ tanh ^pE^ 
sin 29^ tanh ^^E^ + sin 26^^ tanh ^pEj: 



and 



B 



~ AT ^ ' 

k 



-ik.T 



n(ujZ) + niul 



for the interplane parameters 

^^,^+r = E [^os 2^+ tauh - COS 29^ tauh \pE-k 



(13) 
(14) 

(15) 
(16) 
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A' 



i,i+T 



^ ^ [sin 26+ tanh l(3E+ - sin 26]; tanh l(3E, 



(17) 



and 



n[uoi - n[uj. 



and for the carrier number, which must be equal to the number of bosons, and by 
conservation must add to the number of fermions to give unity. 



k 

from the boson part and 

S = ^J2 [cos 2^+ tanh + cos 2^^7 tanh i/JE^ 



(19) 



(20) 



from the fermion part. In (|T3| - ^0|), the sine and cosine factors are given by (^) and 
n{ujk) is the Bose occupation function. 

In this type of model, the bosonic holons occupy only a very small region of reciprocal 
space close to k = 0, and in all respects are effectively condensed at all reasonable 
temperatures. Thus to a very good approximation one may take B = 6 in (|T3|) and 
B' = 5sgn(x') 



m 



so that the bosonic equations need not be considered further. 
It is this approximation which allows us to consider only nearest-neighbour spin order 
parameters in a model with only nearest-neighbour superexchange interactions. Thus 
the problem reduces to a five-parameter one in the fermionic degrees of freedom, with 
the mean-field equations (0), (0), (^), (^) and 

Following Tanamoto et al. |jTT|, |l^, we concentrate on a doping level close to the 
value where the system is close to reproducing the properties of optimally-doped ma- 
terials, and choose 5 = 0.2. Here we expect a solution with finite intra- and interplane 
uniform order parameters over a considerable range of temperatures, and this is borne 
out by calculations. In the low-temperature regime we wish to investigate the nature of 
the transition to a state with singlet order, and the possible coexistence or competition 
of the intra- and interplane ordered states which may exist here. The idea of a singlet 
state with coupled order parameters has been studied in one version of the current model 
by Ubbens and Lee [|T6| , and proposed by Kuboki and Lee [^] as a candidate which may 
explain recent observations of the anisotropic gap in BSCCO materials by photoemis- 
sion [^. A very recent numerical study has investigated the same issue in small. 



bilayer clusters, finding a sharp crossover from predominantly in-plane to predominantly 
interplane singlet order at a value of the interlayer superexchange parameter J± = xJ, 
where x is of order unity, but may under certain circumstances be quite small. We will 
compare our results with those of Ref. |2^ in more detail below. 

In a single-layer system, the gap parameter A may have an arbitrary phase 6, from 
which no physical consequences arise. In a bilayer, the phases of the parameters on each 
layer may differ, thus appearing as A^e*^^ and A^e*^^. When the layers are coupled, there 
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may in general be an additional term in the Hamiltonian of the form H = —e 005(^1 — ^2), 
where the energy scale of the phase coupling parameter is e ~ t^^/W, in which W is the 
bandwidth of the in-plane dispersion. This term may couple to an external magnetic 
field, and gives rise to a variety of Josephson coupling phenomena, some of which are 
integral to the qualitatively new predictions of the interlayer pair tunnelling theory [^|. 



In this work, we assume that the minimum energy state is maintained, so that the phases 
in the two layers remain equal, and both gap functions can be taken simply to be A^, 
as given below (|^). 

In order to make the following study more directly applicable to the physical systems, 
we appeal to the results of Andersen et al. for some additional detail about the 
interplane hopping term. From detailed bandstructure studies on the YBCO system, 
these authors estimate the magnitude to be t\ ~ O.lSt, where in our choice of units the 
maximal band splitting will be %t\. The k dependence resulting from the combination 
of the interplane term with the extended in-plane terms is contained consistently in the 
extended interplane hopping parameters t_L, t'^ and t'| (Fig. 1), and is such that t±(k) is 
small along the FM line in the Brillouin zone. This has been parameterised in the form 

t_L (k) = t\ (cos kx — cos /cy)^ , (21) 

which corresponds to the special case where the bands are degenerate along FM, and 
results from a real-space distribution of the interplane hopping integrals given by t_L = 



32|, = —\t\ and = jt\. The interlayer spin exchange is taken to be nearest- 



neighbour only, and to have a small value J±_ ~ O.IJ, in line with the majority of 
experimental pi, 531 |3H] and theoretical [p2| indications. We note, however, that recent 



infrared transmission and refiectometry measurements fSg] suggest a considerably larger 
value of this parameter in the insulating YBCO system, and so we do not discount this 
possibility. 

To investigate the nature of the spin-ordered state, we consider the linearised gap 
equations (14) and (|17D close to the their transition temperature. In their most general 



form these contain the 3 equations 




+ 



A, 




A. 


cos kx + Ay cos ky 


+ 












A, 


cos kx + Ay cos ky 




|J±A^ 




2C 






A, 


cos kx + Ay cos ky 


+ 


iJ±A^ 




2^ 






A, 


cos kx + Ay cos ky 




|J±A^ 



+ ^ ^ " ' - - tanh 1 (22) 

3 ^-^ / cos kx + cos ky + ^ J_l A_l , i ^ , 





24+ 



tanhi/54 
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Aa; COS + Ay COS ky — ^J±A_ 



tanh ^PC,, 



where = ± \Mk\ and A_l denotes A'. This set of equations can be cast in the 
schematic matrix form 



/ If ^ + I^^ - 1 + L 



jxy I jxy 
-'1 "T -'2 
jx 
-'2 



V If 




Tvy _i_ tw 




(23) 



where = = (cos/cj)^ and J*^ = {cos ki cos kj) a, in which the angled brackets 

( )a denote integration over the Brillouin zone of the energy denominator and thermal 
function of the energy branch a. A coupled transition of all 3 order parameters may 
occur where the determinant of this matrix is zero. The determinantal equation may be 
represented simply by 



a b c 
b a c 
c c d 



(a-b) {a + b)d-2c^ 



(24) 



The in-plane (i-symmetric solution is sought by writing the eigenvector (A^,, Aj^, A_|_) 
as (A, —A, A_|_), where A is real, and it is clear that the only possible solution is A_|_ = 0. 
This is the root a — 6 = of (^4|) above. Thus in the mean-field framework, ci-wave 
singlet pairing in the plane is thus found to exclude any interplane singlet order at the 
same onset temperature T = Tryb- This is the most important qualitative result of 
the current analysis, as it has profound implications for the nature of the singlet spin- 
ordered state which may be realised in this type of model. The in-plane s-wave solution 
(A, A, Aj_) can be found to exist with ^ = — ^, which is the root [ ] = in ( p^ 
above. This is a coupled solution with in-plane s-wave symmetry, which for small J±_ 
admits some admixed interplane singlet order in an amount proportional to the band 
splitting, and at large J±_ represents a predominantly interplane singlet-ordered state. 
We have not been able to find a solution to the system of mean-field equations for the 
in-plane s-wave ordered state for any values of the doping 5 or interplane superexchange 
interaction J±_ in the regime of physical interest, and so do not consider this state further 
here. We will have some further comment on the interplane ordered state below, in the 
context of its possible coexistence with in-plane ci-wave singlet pairing at temperatures 
away from Tjivb- 

Returning to consider the case of a ci-symmetric in-plane singlet order parameter in 
more detail, we note first that if the planes are decoupled (t^ = 0), = and the 
equations for (A^,, Aj^) and A_l decouple. Thus there would be two separate transitions 
for the total order parameter if interplane singlet ordering were not excluded as in the 
coupled case above. To investigate the possibility of coexisting order parameters in some 
temperature regime of the system, one may consider the in-plane d-symmetric ordered 
state and seek a second transition to finite interplane order at a lower temperature. 
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In the absence of interplane singlet order, this state is described by the reduced 
system of 4 equations ( (pTsD , (0), (|16D and (|20|) ), which can be solved readily to yield 
a self-consistent set of solutions for the order parameters as functions of temperature. 
In Fig. 2 we show the transition temperature Tfjy^, representing the onset of in-plane 
singlet order in such a system, as a function of the interplane hopping parameter t^. 
Trvb is little affected by small interplane hopping, but loses most of its value as t° ( pT]) 
approaches J. These results were obtained with a fixed, small interplane spin interaction 

= O.IJ at all values of t° . 

Returning to the system of 5 equations, this can be studied below T^ys at several 
values of (Fig. 2) to seek a transition to non-zero A_l (p!7D. We find that this possibility 
exists only for values of 0.8J < J± < 1.2J: for smaller values the equation cannot be 
satisfied because the coefficient J± of the left-hand side is too small to compete with 
the band splitting, and for larger ones it is likely that the assumption of the interplane 



singlet state being subsidiary to the in- plane state is no longer valid p2[. The range 
of coexistence is quite insensitive to the value of t° , a feature which we believe is 
due to the line nodes in t±k (^)- This result is in agreement with a recent study by 
Eder et al. [^], who investigated spin correlation functions in 16- and 20-site bilayer 
clusters by exact diagonahsation. These authors find for a doped system a very abrupt 
crossover, from a state with predominantly in-plane d-wave spin correlation to one of 
"standing singlets" formed between nearest-neighbour spins in each layer, as a function 
of increasing interplane superexchange J±. The crossover occurs at a value J±_ = xJ, in 
which X varies over a parameter-dependent range of values below unity. That interplane 
singlet order is not compatible with in-plane (i-symmetric singlet pairing was also found 
by Liechtenstein et al. in a bilayer antiferromagnetic spin fluctuation model, which 
was shown to favour the existence of the anisotropic s-wave state. 

These results raise the interesting question of how the (i-symmetric in-plane singlet 
order in an isolated plane may evolve as interplane hopping and superexchange increase. 
One alternative to the admixture of real order parameters considered above, which has 
been shown not to evolve below some significant threshold value of the coupling, is the 
development of imaginary components of the order parameter. If the order parameter 
has the form A^^ = Ai + iA2, A^, = — Ai +iA2, one may formulate a system of equations 
analogous to ( ^2]) above, in which the real part is the single-layer equation for in-plane 
(i-order and the imaginary part a coupled equation for A2 and a purely imaginary 
interplane component iA±. This has been termed the d + is state, and would give a 
completely nodeless gap because the order parameter has a finite imaginary part where 
the real part vanishes. However, the real and imaginary parts are decoupled, so there 
is no reason to expect the two transitions to occur simultaneously. A second transition 
below TuvB to a nodeless state has been proposed as the explanation for a feature in 
the NMR spin-echo relaxation rate [^, and so this state is of both experimental 
and theoretical interest. We have sought a transition to this state in the same way as 
described above, and again found that it did not exist within the present framework 
below values of Ji ~ J. 
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A different proposition for the existence of the coupled singlet state is given by the 
argument that fluctuations of the gauge- field in a bilayer system are considerably 

more effective in suppressing in-plane d-wave order than interplane s-wave order, and 
thus act to allow the latter to survive in the presence of the former. The same authors 
also argue for enhancement of the interlayer coupling J±(q), to values which may fall 
in the coupled regime, close to the antiferromagnetic wavevector {tt,tt). Finally, if the 
system is orthorhombic, i.e. 7^ ty and Jx 7^ Jy in the plane, one may show ^ 



that the possibility of a coupled transition does exist, and thus that a mixed "rf + s" 
or "(i + is" state of coexisting intra- and interplane order parameters is justifiable in 
the physical parameter regime. However, there is as yet no definitive evidence that the 
BSCCO system, whose photoemission spectra have been interpreted to suggest a d + s 
state [pri, shows any orthorhombicity of a type which would couple appropriately to 



both gap parameters. 

The range of J± for which a mixed state of in- and interplane singlet spin order may 
exist is unrealistically large for the high-Tc superconducting materials, and in addition 
is unphysical in the context of a model where the superexchange interaction should be 
the result of a second-order hopping process, with magnitude J± ~ In order to 
investigate the properties to be expected from the model, and to compare them with 
experiment, we wish to restrict the discussion to the physical parameter regime, which 
will require considerably smaller values of J±. We will thus proceed by taking the 
interplane singlet order parameter to be zero, and to study the 4-parameter problem 
for the uniform in- and interplane order parameters x ^^^^d x'j the in-plane spin gap 
parameter A and the chemical potential fi, albeit in the presence of a finite interplane 
spin coupling constant. 

Following Andersen et al. we now fix the magnitude of the interplane hopping 
term at t° = 0.15t. We choose to set the interplane spin coupling to the small value 
J± = O.IJ, noting that if in fact this parameter is found by more detailed experiments 
to have a somewhat larger value, its role in the structure of the current theory is such 
that it has a negligible effect on the mean-field solutions and Fermi surface shape. In the 
spirit of the single-layer analysis of Tanamoto and coworkers |]TT|, ITTl , we wish to choose 



the extended transfer integrals t' and t" in such a way that the Fermi surfaces of the 
bonding and antibonding bands reproduce as many features as possible of the physical 
bilayer systems of interest, BSCCO and YBCO. As was shown clearly by these authors, 
the magnetic properties of the system are strongly influenced by the exact shape of the 
Fermi surface. 

As a prototype of the BSCCO system, in Fig. 3(a) we show Fermi surfaces in the 
first quadrant of reciprocal space for the parameter choice t = 4 J, t' = — |t and t" = |t. 
For the fixed value of this leads to a pair of dispersion surfaces with fiat regions 
around the (vr, 0) saddle points, the important feature of which is that the "extended 
saddle point" regions of the bonding band lie below the chemical potential, while those 
of the antibonding band lie above it. Thus the bonding (outer) band has the open 
shape characteristic of YBCO-like systems |T^, while the antibonding (inner) band has 
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the closed, LSCO-like shape; we will use the terms "open" and "closed" to refer to the 
nature of the Fermi surface in the first Brillouin zone, in an electron picture. In addition, 
the parameters are chosen so that the antibonding band is rather strongly nested, by 
which is meant that the Fermi surface has significant regions lying almost on a straight 
line parallel to g^; + = 0, so that the same wavevector q can span many pairs of 
Fermi surface points; in detail, this surface is slightly concave in relation to the given 
line. Finally, the bands appear to be degenerate along the direction FM in the Brillouin 
zone, in accordance with the calculations of Ref. [^, and so we retain the form for t±k 
(|2l|) proposed by these authors. These are the features of the published interpretation 



of the photoemission data from BSCCO in Ref. and so we take this parameter set 
to characterise the BSCCO system in the remainder of this study. 

To describe the YBCO system, in Fig. 3(b) we show Fermi surfaces in the first 
quadrant for the parameter choice t = 4 J, t' = —^t and t" = jt. In this case the 
parameters are such that the extended saddle point regions of both bands are below the 
chemical potential, so that both have open Fermi surfaces. However, the antibonding 
band is very flat in these regions, and the Fermi surface shows a rather narrow "neck" , so 
is close to becoming closed on small changes of the parameters. We emphasise in passing 
that these strong alterations in Fermi surface shape are a consequence of very small 
changes in the parameters t' and t" , because of the fiat nature of the saddle-point region. 
Here we require that the bands are not fully degenerate in the FM direction, and so 



alter the form of the interplane term to t±k = | {cosk^ — coskyf + 2 . This alteration 
corresponds to the case where the interplane hopping integrals are given hy t±_ = 0, t'j^ = 
andt'_[ = or in other words that the interplane hopping is less extended in this 
case. The function in reciprocal space remains minimal along the FM. The fitted features 
are those of the YBCO system, as calculated from first principles in Ref. and are 
fully consistent with results from the most accurate photoemission studies of YBCO 124 
PB| and 123 to date. Henceforth we take this parameter set to give the prototypical 
YBCO system for analysis of the physical properties. We believe that the relatively 
large value of t" = is reasonable, as it remains consistent with the expectation from 
microscopic models for the transfer of extended quasiparticles in the Cu02 plane [p9| . 
The Fermi surfaces shown in Fig. 3 are independent of the perpendicular wavevector kz, 
which does not appear in the quasiparticle dispersion (||) due to the assumption of no 
coupling between unit cells. In the slave-boson decomposition at the mean-field level, 
because the holons are concentrated at g = 0, calculated photoemission spectra are 
exactly characteristic of the spinon Fermi surfaces (Fig. 3), and contain no alterations 
reflecting any effects of the charge degrees of freedom on the spin. 

With the above choices one may solve the mean-field equations of the 4-parameter 
fermion system, to obtain all of the parameters as functions of temperature. In Fig. 4 
we show the form of the in-plane singlet spin order parameter A(T) for both BSCCO- 
like (Fig. 4(a)) and YBCO-like (Fig. 4(b)) systems. A shows a BCS-like second-order 
transition, and the transition temperature T^ys sets the characteristic scale for spin- 
related properties in this type of model. We turn now to the computation of two such 
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properties, the dynamic susceptibility and the phonon anomahes. 



3 Dynamic Susceptibility 

The dynamic susceptibihty x(q, g^jCj) contains all information about the spin response 
of a magnetic material, at all frequencies and wavevectors. Its zero-frequency limit is 
sampled at certain wavevectors by NMR experiments, while the full dynamical quantity 
is measured by inelastic neutron scattering, as discussed in more detail from the view- 
point of this type of model in Ref . |]T^ . Here we will compute the dynamic susceptibility 
in the random-phase approximation (RPA), one method by which spin fluctuation en- 
hancement may be taken into account, and illustrate the results as scans of in-plane 
wavevector q or frequency uj with the other variable held fixed. We will show also the 
variation of the susceptibility with the wavevector component g^, and explain the results 
in terms of intra- and interband scattering processes. 

In a bilayer system, the imaginary part of the bare, retarded susceptibility XoIQ) 1zi ^) 
is a 2 X 2 matrix which can be calculated by analytic continuation of the associated 
thermal function 

, . , / (S^.S^) (S^.S^) \ 

in which (S^S"^) denotes (S'(q, g^, ico'„).S™(q, g^, zti;„)). The spin operator S may be 
written in terms of the operators /fc, transformed to the quasiparticle operators 7^ and 
the resulting expressions evaluated in the diagonal basis. Omitting details of this lengthy 
procedure, the bare susceptibility assumes the form 

Xo(q,g.,^)-(^^-.,..^_^q^^) x4q,-) j' ^^^^ 

where the diagonal and off-diagonal parts are quantities whose imaginary parts have the 
form 

X'4(q, 0;) = E E [CtJt, - El + Et^^) -5(uj + El- E^^^ 

a/3 ° k 



+ C,-,T,-^ [5 (u + El + -5(u:-El- E^^^)]] ,(27) 



and similarly for the real parts x± ||T3l- Here the indices a and [3 denote the bonding 
or antibonding bands, = 1 and A~ = 26af3 — 1, 

Cj^ = l±cos2(0,^^-^f) (28) 
gives the coherence factors and 

= /(±Ef ) - f{E^^^) (29) 
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gives the thermal factors. Thus each component of xo is a combination of 16 terms, 8 
for normal scattering processes and 8 for anomalous ones. All of the contributions are 
additive in x+, but the interband contributions (a 7^ P) are subtractive in x_. 

Note that the only dependence on is that contained explicitly in the phase factors 
of the off-diagonal terms in (p6|) . This phase will appear in the imaginary part as the 
legitimate coefficient of terms proportional to a 5-function in energy which are respon- 
sible for dissipative processes, and is not forbidden by the symmetry of x"- On taking 
the symmetric combinations required to compute any physical quantity, the phases will 
combine to give a real result. 

The RPA susceptibility is a sum of chains of bare polarisation parts connected by 
the spin interaction J(q, g^), and can be written in the form of a matrix Dyson equation 
as 

= xo(g,^) - xo(g,^)J(g)x(?,^) (30) 
= (i + xo(g,^)J(g))"^xo(g,^), (31) 

where q denotes (q, g^). The matrix form of the spin interaction on the bilayer is 

3(aq)-( ^11^^) e^^^^^^^x(q) \ .33^ 
Jlq,g.J - J||(q) j' ^^^^ 

in which by the assumption of nearest-neighbour interactions only, one has J||(q) = 
J{cosqx + cosgy) and J±(q) = J±. 

Omitting the functional dependences and superscripts for clarity, the RPA suscepti- 
bility from (^) has the form 



^ \det\ [ e-^'^ {x- + Mxl - X'-)) X+ + Mxl - x'^) 

where denotes g^rc and the determinant of the matrix {1 + xoJ) has the separable 
form 



\det\ = l + 2[x+Jii+X-J±) + [xl-X-)[J^-Jl 

= {l + iJi\ + J±)ix+ + X-)){l + iJ\\-J±)ix+-X-))- (34) 

The measured quantity in an experiment such as the determination of the neutron 
scattering cross section is the sum of all components of the imaginary part of the sus- 
ceptibility matrix, i.e. 

S'M'^Ef^, (35) 

where and u denote the layer index. Thus the diagonal terms in (^) are the effective 
RPA susceptibility for spins in a single plane of the bilayer, and the off diagonal terms 
give the interplane contribution, which will be proportional to cos g^r^. 
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By performing the summation and separation, the final form of the cross section is 

^ J \ (X+ + x^)cos^|9.^c , (x+ - X-) sin^ |g,r, | 
b"\q, uj) (xlm { — — - — ^ + — — - — ^ ) , (36) 



a more general version of the result given by Ref. ||T6|. Returning to the equation 
(|27|), one observes that the components may be cast in the form x± = Xaa ± Xai3 ± 
Xi3a + where each of the terms Xaa' has contributions from 2 normal and 2 anoma- 
lous scattering processes, and is either purely intra- {a = a') or interband. The 2 
terms in ( P^ ) are then either purely intraband (x+ + X- = '^{Xaa + X/?/?)) interband 
(x+ ~ X- = '^{.Xap + X/3a)) ill nature. Thus the wavevector may be used to select 
which type of scattering process is measured, so that the contributions of each may be 
compared. 

In Figs. 5-7 we give results for the wavevector and frequency dependence of the 
imaginary part of the dynamic susceptibility x ~ "Yli^u X/iiy 

for the BSCCO- and YBCO- 

like systems, based on the Fermi surface parameters chosen in section 2. Here we will 
show only results for the system at low temperatures (specifically, T = O.ITjivb), so that 
it is in the singlet RVB state; the qualitative features of the evolution from the high- 
temperature state are well illustrated in Ref. |T^. In Fig. 5(a) is shown x"(ct; ^) for the 
BSCCO system as q is scanned along the symmetry directions of the Brillouin zone. The 
frequency is fixed at = 0.3 J, above the value of the twice the maximal gap A^, and the 
out-of-plane wavevector at = so that only scattering processes within the bonding 
and antibonding bands may contribute. One observes a significant incommensuration 
due to the strongly-nested shape of the antibonding band, and a general increase in the 
scattered intensity for values of q near Q = {tt,7c), the antiferromagnetic wavevector 
which spans many points close to the Fermi surface of the bonding band (Fig. 3). This is 
a combination of the Fermi-surface features documented in Ref. [|17| for both closed and 
open types, which clearly vary little with energy transfer. In Fig. 6(a) is shown the same 
quantity for = vr (in "bilayer units", Vc = 1), so that only bonding to antibonding band 
scattering processes are sampled. Here the scattered intensity is considerably higher, 
indicating a greater enhancement from the relevant denominator in (|36D, and is largely 
commensurate as a result of poor nesting between the two bands. In Fig. 7(a) is shown 
x"(q, cij) at fixed wavevector q = Q as the frequency is scanned up to values of order 
J, for 3 choices of Qz and at low temperature. The results for interband scattering are 
very similar to those obtained in Ref. |]T^ for the single layer with open Fermi surface, 
in that there is a strong peak due to states pushed out of a gap (Aa(ka) + A^(k^) here, 
with ka — k/3 = (vr, tt)), followed by a broad region of spectral weight up to a shoulder 
at u; ~ J which characterises the maximal energy separation of significant numbers of 
points whose wavevector separation is (vr, vr) for the energy dispersions of the system. 
For intraband scattering, the shoulder appears at a similar energy, but the peak feature 
is missing; this can be understood due to the fact that there is no contribution from 
processes involving the antibonding band, all of whose spanning wavevectors are smaller 
than Q, with the result that the denominator in the first term of (|36|) is far from being 
small, and so provides very little enhancement. 
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Fig. 5(b) shows the variation of x" with q in the Brillouin zone at fixed frequency, low 
temperature and Qz = for the YBCO system, where it is clear that the spin excitations 
are commensurate, and strong over a broad range of q around (vr, vr). This result appears 
from Ijl^ to be characteristic of the open Fermi surface shape; the only qualitative effect 
of the higher energy transfer shown here in the singlet RVB state is that one no longer 
observes the peaks which arise due to scattering processes between the gap nodes. For 
= TT (Fig. 6(b)), the peak is again commensurate, as is to be expected as the 2 bands 
involved both have the open shape, and once again has higher intensity than that from 
intraband processes. Fig. 7(b) shows that the excitation spectrum as a function of u 
once more has a shape similar to that seen in the single-layer case, with the exception 
of clear evidence of 2 separate contributions to the peaks at each q^. This is the result 
of having 2 pairs of Fermi-surface points separated by ('/r,7r), so that the characteristic 
gap combinations in each case may differ, and is most pronounced in the = vr case. 
The shoulder feature for the YBCO bands occurs at o; ~ 1.2J. 

We may conclude by noting that at the RPA level the dynamic susceptibility is very 
much dominated by Fermi-surface scattering processes, which are strongly favoured over 
even those involving nearby points, and in this sense the results show little qualitative 
difference from those of the single layer. The calculated result that interplane contri- 
butions are stronger than intraplane ones over most of the energy range (Fig. 7) arises 
primarily from their smaller enhancing denominator (^6]), and is consistent with the 
results of experiment in YBCO [^. In the absence of coherent interplane hopping [p5| , 
a different explanation of this modulation would be required. However, with regard to 
overall agreement with experimental observation, we do not consider the results given in 
Figs. 5-7(b) to be a satisfactory description of the measured magnetic response of YBCO 



materials Recent experiments show few features in the spin excitation spec- 

trum other than a strong and remarkably narrow peak at 41meV in the superconducting 
state, although there is evidence for some broad spectral weight below this energy in 
the normal state. The peak occurs only close to the wavevector Q, and does not seem 
to disperse significantly with q or u, so that the band picture given above, and pro- 
posed by many authors, would appear to be a poor candidate for a full explanation. 
This is in contrast to the spin spectrum of the single-layer LSCO material ^ 
where the current form of theory gives a good account of the experimental features. In 
YBCO, the narrow and non-dispersive nature of the 41meV peak suggests the existence 
of a resonance, and there has been a theoretical proposal for this within the class of 
tight-binding, strongly-correlated systems in two dimensions |^^. A full investigation of 
this possibility in the present framework will be the subject of a subsequent publication. 
Finally, we comment that no measurements of dynamic susceptibility have so far been 
performed on BSCCO materials, because these have proven extremely difficult to grow 
as high-quality single crystals of the required size 
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4 Phonon Anomalies 



The coupling between the spin and lattice degrees of freedom in the t-J model has been 
investigated recently in Ref. [l^. These authors use the modulation of the parameters tij 
and Jij by the c-axis oscillation of in-plane O atoms in a single-layer model to explain the 
superconductive phonon anomalies in the important and well-characterised 5-symmetric 
modes of the YBCO plane, as well as to obtain qualitative agreement with the results 
for other materials classes, with the observation of the spin gap in underdoped YBCO 
compounds, and with isotope effects at optimal doping. The reader is referred to this 
paper for the background to this section. 

In Fig. 8(a) are shown the four c-axis modes of planar O atoms in the bilayer system: 
at the single-layer level there is no difference between the 340cm~^ mode and 
the 193cm~^ mode, or between the 440cm~^ Aig mode and the 307cm~^ mode. 
While this may be realistic for the i?-symmetric modes, which develop no net dipole 
moment, it was found not to afford an acceptable description of either the mode, 
where there will be significant interplane charge motion during the phonon oscillation, or 
the Aig mode, in which charge may be transferred out of the bilayer unit to neighbouring 
planes. Here we investigate the phonon anomalies in the bilayer in the presence of finite 
interplane coupling and superexchange, which as above lead to a finite interplane uniform 
spin order parameter, but no interplane singlet order. There have been indications 



within a phenomenological weak electron-phonon coupling model |^ that interlayer 
hopping terms are sufficient to reproduce the phonon anomalies measured in the 
infrared- active A2U mode. 

^ , in the bilayer it is necessary to recast the phonon-spinon coupling 



Following Ref. 



vertices in terms of the diagonalising quasiparticle operators. Each operator is a sum 
of 4 operators with coefficients given by (^), so that each vertex is a sum of 16 terms, each 
involving a contribution from both spin orientations. For the "uniform RVB" vertex one 
obtains on substitution 



.l(2)t .1(2) 



CiTifcTifc + Si-rlklLk T Ci27ifc72fc T SulUll-k 

+S'i7i_fe7ifc + Ci-fi^kll-k =F 5'i27i-fc72fc ± Cuji^kll^k (37) 
TCi2llklik T SullklLk + C'272fc72fe + '^272fe7Lfc 
=F5'i272-fc7ifc ± C'i272-fc7i-fe + 5'272-fc72fe - C'272-fe7Lfc ^ 



and for the "singlet-RVB" vertex 



-5171^71^ + Ci7ifc7i_A: ± 51271^7^ T Ci2jikl2-k 
+Ci7Lfc7ifc + 5i7j_fc7i_fe T Cull^kllk T 5i27Lfc72-fe(38) 

±5'l272fc7L =F C'l272fc7l-fe - S2j2kllk + C2l2kl2~k 

T Ci2il_ki\k T Si2il_kli-k + C2il_kilk + '527Lfc72- J , 
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where Ca and Sa denote cos 29'^ and sin 2^", and C12 and S\2 are the combinations C1C2 — 
S1S2 and C1S2 + S1C2, where Cq, and denote cos^^ and sin 6'^. The complexity contained 
in 2x2 matrix vertices and Nambu propagators for the fk problem is transferred solely 
to the 4x4 matrix vertex in the quasiparticle problem, conserving the number of degrees 
of freedom. 

In the even modes {Aig and -Big), the motions of the O atoms are out of phase 
between the planes, so that the effect of their motion on t and J is the same in each 
plane. Defining the phonon coordinate as positive when the oscillating O atom moves 
towards the plane of the Cu atoms (Fig. 8(b)), one has in obvious notation u]. = Uy = 

'^l = '^l = for ^19 '^l = = = = for ^ig- Thus in both cases 
the effective vertex from summing the contributions of both planes in ( ^7|) and ( ^8l) 
above is to cancel all of the terms mixing the indices 1 and 2 for both x and y bonds, 
leaving only those terms related to intraband processes. Conversely, the odd modes 
{ul = Uy = —ul = —Uy = u) and i?2« [u], = —Uy = —ul = Uy = u) will have 
vertices coupling only to interband scattering processes. This very simple but profound 
separation of the contributions may be understood readily, in that the even modes 
produce a pattern of modulations of t and J which is even in the unit cell, leading only 
to intraband transitions, and conversely for the pattern of modulations arising from 
the odd modes. The lowest-order quasiparticle polarisation terms contributing to the 
phonon self-energy may then be calculated to deduce the frequency shift and linewidth 
broadening of each type of mode. For the even modes (intraband), the contribution 



from normal scattering processes will vanish in the limit g — > 0, and one may follow |jT9| 
exactly to obtain 



tanh i ^] 1 tanh ( ^ 



,2^y p+_ 1 I 2T ; „_ 1 ''''''' \ 2T 



SuJ = c{Xjjy-Y.Fjt -—-2 ^^ + Fk , 72 (39) 



' is 



Here c = (£) (u^) = 1.18 x 10-^ (u^) = = (O.OSSA)^, in which uo = 340cm- 

taken from the Big phonon frequency and M is the mass of the O atom, and E"^ are 
given by (0). The form factors are 



Ft = 



Big 



F^ = 2ASt i^k±\Mk\) + 'ixik' Aig, (40) 



where 'fk = cosk^ + cosky, rjk = cosk^ — cosky, A = (Ajj) and x = iXij) + 37- I^i 
this calculation we have not included an "interplane uniform" vertex with contributions 
proportional to B[ and x^, and will discuss the effects of such terms in detail below. The 
linewidth broadening can be computed from an analogous expression for the imaginary 
part of the self-energy. 

For the odd modes, the interband scattering expression and form factors become 
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more complicated. One obtains 

5uj 



El. 



o 1 2fE, 



2 Et + El 



{Et + El 

in which the form factor for normal scattering events is 



2 (/(^.") - f^El 



(41) 



y2A7fc(ciS2 + S1C2) + XVk{ciC2 - S1S2) 
V2Ar]k{ciS2 + S1C2) + X7fc(ciC2 - S1S2) 
and that for anomalous scattering events is 



Fk 



F^ 



A/2A7fc(ciC2 - S1S2) + XVk{ciS2 + S1C2] 

V2Ar]k{ciC2 - S1S2) + X7fc(ciS2 + S1C2) 



B2U 

^2n; 

B2U 
A2u- 



(42) 



(43) 



It is straightforward to show that in the limit where the bands become degenerate, both 



sets of expressions reduce to the result expected from |]T9| for the single layer. 

In Figs. 9-12 we show the phonon frequency shift 600 calculated using the above 
formulae, and the linewidth broadening ST, which is calculated from analogous expres- 
sions. We compute these anomalies only for the YBCO system, where the buckling 
of the O atoms in the bilayer structural unit is uniformly directed along the c-axis to- 
wards the symmetry plane of the unit cell, an arrangement which may yield constructive 
combination of the vertex contributions to give a linear spin-phonon coupling term in 
cases where the symmetries of the phonon distortion and the superconducting gap are 
compatible. By contrast, in BSCCO the single layer is orthorhombically buckled, into 



a tilting distortion which is odd in the unit cell [47|, so there is no linear coupling term 



and the anomalies are expected to be negligibly small [|l^], as confirmed by experiment 



In Fig. 9(a) is shown 6uj as a function of frequency for the (/-symmetric modes, at 
g = and at a temperature well below the singlet ordering transition Trvb (Fig- 4). 
As in the single-layer case, low-frequency 5-symmetric phonons will be softened by the 
spinon contribution, while high-frequency modes will show hardening. The character- 
istic frequency of the sign change is set by the value of the gap near the heavily- 
contributing (vr, 0) points for the intraband processes in question. Quantitatively, it is 
somewhat smaller than in the single-layer solution, requiring that the analogue of the 
maximally-shifted Big mode in this system appear at ~ 0.2J, in comparison with the 
experimental mode frequency cuq — 0.3J. The imaginary part (Fig. 10(a)) has the ex- 
pected form of a peak at the frequency where the real part changes sign, and shows that 
only phonon modes with frequencies close to this value will be significantly broadened. 
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For even, 74-symmetric modes one finds also the result of the single-layer system, that 
the anomalies are negligible. The symmetry effects which cause the linear contribution 
to be small in this case are not affected by the introduction of interplane terms. 

For odd (w-symmetric) modes one sees a different feature (Figs. 9(b), 10(b)). Because 
these couple only to interplane scattering processes, the characteristic frequency which 
determines the change from phonon softening to hardening is set not by twice the gap 
Afc, but by twice the interband energy separation ~ 6t_i_k, again at the (vr, 0) points 
where it is largest. From the results of Ref. which were used in setting the interplane 
hopping parameter and are consistent with available experimental evidence, this 
separation is of order J, and so appears on an altogether different energy scale. Thus 
all M-symmetric phonons in the usual energy range will soften at the singlet ordering 
transition within such a model, and the linewidth broadening will be negligible. 

In Figs. 11(a) and 12(a) are shown the temperature dependences of 6uj and for the 



5f-symmetric modes; this is the quantity measured by experiment I^T], and should 



be directly comparable with the results for the Big and Aig modes. One sees again that 
the results are very similar to the single-layer case: for the Big mode there is a gradual 
development of the frequency shift over a range of T below its onset at Trvb, and a peak 
in the linewidth broadening there, both in good agreement with experiment ||5l]]. The 



chosen mode frequency here is luq = 0.2J, which is somewhat below the experimental 
value, but is consistent with the location of the Big phonon being that for a maximal 
anomalous effect, and remains in a reasonable degree of correspondence with the physical 
system given the nature of the model. For the Aig mode, the chosen frequency should 
be in the phonon hardening regime, but we show again uo = 0.2J to emphasise that 
there is very little effect, and we comment on this in detail below. In Fig. 11(b) is the 
T dependence of 6uj for a mode of frequency uq = 0.2J and B2u symmetry. Here uq is 
so far from the characteristic frequency of the changeover from softening to hardening 
of M-symmetric phonons (Fig. 9(b)) that the development of the anomalous frequency 
shift is effectively immediate at T^ysi as seen less obviously in the single layer model. 



and in good agreement with neutron scattering studies of this mode in YBa2Cu307 [p3 
Similarly, ^F for such phonons as a result of coupling to the spin is negligible within 
this model. 

We have considered only g = phonons here, for comparison with the results of 
Raman and infrared spectroscopy. For finite in-plane wavevector q, one may expect the 
results to be similar to the single-layer case [l^ for both g- and u- symmetric modes, 
namely that the anomalies show some features of the Fermi surface shape, and fall to 
zero as q approaches the (vr, 0) and (vr, vr) points. However, at finite the intra- and 
interband scattering processes become mixed, so that both mode symmetries may show 
both characteristic energy scales. Of particular note is the situation at q = 0, = vr, 
where the anomalies will remain large and the contributing processes will be completely 
reversed in (intra- or interband) nature. Thus it should be possible by neutron scattering 
to see that the onset of the frequency shift 5uj for the Big phonon is very much sharper 
in temperature, if smaller in overall magnitude, as its frequency will no longer be close 
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to the relevant energy scale. Similarly, no alterations in the linewidth of this phonon 
should be detectable. 

We have shown that the bilayer system retains a good description of the main qual- 
itative features of the 5-symmetric modes of planar O in the bilayer structural unit, 
although we stress again that the g- and u-symmetric modes have very different char- 
acteristic energy scales in the coupled bilayer system. However, the current level of 
treatment does not contain the physics required to explain the experimental results for 
the y4-symmetric modes. Here we have included the interlayer hopping only through 
the way in which it arises in the quasiparticle transformation, and in this sense tx^ is 
treated on an equal footing with the terms t'^. and t'^, which are taken to determine 
but not to contribute qualitatively to the phonon anomalies |]l9l. This is also the 



approach taken by Ref. but these authors chose the value of v\_ to coincide with the 
phonon frequency, in order to enhance the anomaly in their approach, and in addition 
used both less restrictive symmetry constraints and an isotropic gap function Aq. 

We have neglected the possibility of a contribution from modulation of the magnitude 
of the parameter t^, because in none of the phonon modes considered do the Cu atoms of 
the bilayer change their relative separation. If one wishes to include such a modulation, 
for example on the grounds that additional transfer may occur through the neighbouring 
O atoms (whose separation changes in (^-symmetric modes), it would appear as a term 
-^±Xfc (/fcVfc + h.c.^ in the Hamiltonian. Substitution of the quasiparticle operators 
yields the same type of vertex terms as in the expressions (^) and but with 

considerably smaller coefficient for reasonable values of the interplane phonon coupling 
Aj^. Thus the additional term would neither make a significant contribution nor appear 
with a symmetry factor which combines constructively with A-symmetric phonons. For 
this last reason, even a very strong modulation of the coefficient f\_ would not explain 
the experimental observations within the current theory. Such a modulation may in fact 
be present in the mode, where a very strong dipole moment is developed during 
the oscillatory motion, but it appears that this must be described by additional physics, 
related to the strong effects of the motion on interlayer charge transfer. Similarly, in the 
Aig mode one must also take account of charge transfer, which may occur not between 
the planes of the bilayer but between each plane and the neighbouring BaO layers 
(in which 0(4) is the apical oxygen atom of planar Cu(2)). Recent Raman-scattering 
experiments suggest that the 440cm~^ Aig planar mode and the 500cm~^ Aig mode 
of apical O behave as split levels of a single, fundamental phonon frequency, indicating 
very strong coupling between in- and out-of-plane modes of these symmetries. 

We conclude this section by noting that apart from the change in the nature of the 
M-symmetric modes discussed in detail above, none of the qualitative features extracted 
from the single- layer model |19| are altered in the bilayer system. In particular, the 
important result that the spin gap in the low-doping regime of the model is observable 



due to spin-phonon coupling |^ is retained, a point which reflects that the current 
model does not contain a difference in such properties between mono- and bilayer oxide 
materials. The isotope effect arising as a result of the coupling between the lattice and 
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the spin sector, which is responsible for superconductivity, is identical to that deduced 
from the single-layer model, because the bilayer treatment also does not reproduce the 
anomalous contributions from A-symmetric modes. 



5 Summary 



In summary, we have analysed the spin properties of a bilayer CUO2 system, using 
the mean- field, slave-boson treatment of the extended t-J model. We find that the 
existence of d-symmetric spin singlet order in the plane, favoured by the single-layer 
model, acts to suppress interplane order completely below a characteristic, finite value 
of the interplane coupling. At this coupling, there is a very rapid crossover to a state 
of pure interplane singlet spin order, suggesting that the two types do not coexist. 
Working within the physical regime of interplane coupling, which mandates exclusively 
in-plane singlet pairing, we show that the observed Fermi surfaces of both BSCCO- and 
YBCO-like bilayer systems may be reproduced with realistic extended in-plane transfer 
integrals. 

For each type of system, we have calculated two microscopic and physically observ- 
able properties related to the spin sector. The dynamic susceptibility x(ct, t^), computed 
at the RPA level, shows the spin fluctuation spectrum at all frequencies and wavevec- 
tors; the commensuration and excitation properties may be explained simply in terms of 
Fermi-surface scattering processes, but the remaining discrepancies between the theory 
and experiment in YBCO reveal that the problem is significantly more interesting and 
complex than mean-field physics. Anomalies in the frequency and linewidth of phonon 
modes in the YBCO system show a strong contrast between even ((^-symmetric) modes, 
which couple to intraband scattering processes, and odd (w-symmetric) modes, which 
couple only to interband scattering. The model provides a reasonable account of the 
S-symmetric modes, which have no net electric dipole moment in each layer and so are 
not accompanied by charge motion, but cannot describe the A-symmetric modes, which 
do not satisfy this condition. 

We conclude by commenting that because in this model the physical parameter 
regime remains dominated by the d-symmetric in-plane singlet order of the single-layer 
case, the results for observable quantities in the bilayer system differ qualitatively from 
those of the single layer in only a small number of features. We have carried out semi- 
quantitative calculations which may be compared in detail with experiment. This makes 
it possible to quantify the strengths and limitations of the current t-J model formulation, 
to learn which properties appear due to the bilayer nature, and to gain some evaluation 
of the restrictions of the mean-field theory and RPA approaches. Thus we may deduce 
the directions in which to look for the additional physics which will be required for a 
more complete description of the high-Tc systems. As an example of this process, we 
have assumed the existence of coherent single-particle transfer between the planes of the 
bilayer, and investigated the experimental consequences. The interlayer pair tunnelling 
theory pSf argues that such transfer cannot be coherent, and leads to some very different 
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physical predictions, which remain to be computed for similarly detailed comparison. 
Continuing rapid improvements in crystal quality, and in a variety of experimental 
techniques, will provide crucial evidence in favour of one of these viewpoints in the near 
future. 
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Figure Captions 

Fig. 1: Schematic representation of the bilayer system, showing transfer integrals for 
quasiparticle hopping and spin superexchange interactions. 

Fig. 2: Variation of singlet spin ordering temperature Tbvb with interplane hopping 
parameter tj^, at fixed J± = O.IJ. The in-plane transfer integrals are t = 4 J, t' = —^t 
and t" = \t. 

Fig. 3: Bonding (o) and antibonding (x) Fermi surfaces for bilayer systems with doping 
5 = 0.2, shown in the first quadrant, (a) BSCCO-like system: parameters are chosen, 
following Ref. to give one open and one closed Fermi surface, strong nesting of the 
latter (antibonding) surface, and band degeneracy in the PM direction, (b) YBCO-like 
system: parameters are chosen, following Refs. P^, EBl to give open Fermi surfaces 



for both bands, but with that for the bonding band close to "closing", and to lift the 
degeneracy in the PM direction. 

Fig. 4: Self-consistent solutions for the singlet-RVB order parameter A(T). (a) BSCCO- 
like system, (b) YBCO-like system. 

Fig. 5: Dynamic susceptibility calculated as a function of q in the 2d Brillouin zone at 
fixed frequency u = 0.3 J, low temperature and Qz = (intraband processes only). In our 
notation, the X point corresponds to (vr, 0) and is symmetrical with Y, while M denotes 
the point (7r,7r). (a) BSCCO-like system: note strong incommensuration arising from 
nesting of the antibonding band, (b) YBCO-like system: note commensurate nature of 
the peak intensity. 

Fig. 6: Dynamic susceptibility calculated as a function of q in the 2d Brillouin zone at 
fixed frequency u = 0.3 J, low temperature and = vr (interband processes only), (a) 
BSCCO-like system: note that intensity of interband processes is considerably higher 
than in the intraband case, and that they exhibit only minor incommensuration. (b) 
YBCO-like system: susceptibility remains commensurate at all g^. 

Fig. 7: Dynamic susceptibility calculated as a function of uj at fixed in-plane wavevector 
q = (vr, vr), with = (x), = tt/2 (□) and = vr (o). (a) BSCCO-like system. 
The main peak is due to states pushed out of the gap, and is absent for the g^ = case, 
primarily because the antibonding Fermi surface is everywhere smaller than (vr, vr)). The 
shoulder at ~ J indicates the energy separation above which there are no more states 
separated by (vr, vr) in q space for the BSCCO-like dispersion, (b) YBCO-like system. 
The two contributions to each peak are due to gaps with different values on each band 
at the wavevector k corresponding to (vr, tt) scattering. For the YBCO dispersion, the 
"shoulder" occurs at a; = 1.2J. 
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Fig. 8: (a) Atomic motions and experimentally observed frequencies of the four phonon 
modes involving c-axis motion of in-plane O atoms in the buckled Cu02 bilayer. (b) 
Phonon coordinate labelling convention. 

Fig. 9: Phonon frequency shift 5ui as a function of frequency and at T = 0.05TRys 
in the YBCO-like system, (a) Big- (o) and Ai^-symmetric (x) modes: the frequency 
of the sign-change is characterised by the gap value at the strongly- contributing (tt, 0) 

points, because only intraband processes arc involved, (b) B2u- (o) and A2u-symmetric 
(x) modes: because the phonon coordinate couples only to interband processes, the 
frequency of the sign-change is characterised by the value of the interplane hopping 
term at (tt, 0) points. 

Fig. 10: Phonon linewidth broadening 5r at g = 0, as a function of frequency and at 
T = O.OSTrv-b in the YBCO-hkc system, (a) Big- (o) and ^ig-symmetric (x) modes, 
(b) B2u- (o) and yl2u-symmetric (x) modes. 

Fig. 11: (a) Phonon frequency shift 5uj as a function of T for hypothetical Big- (o) and 
Aig-symmctric (x) modes of frequency uoq = 0.2 J, at g = for the YBCO-like system. 
Because the mode frequency is close to the characteristic frequency of the sign-change 
(Fig. 9(a)), the full frequency shift develops over a range of T. (b) Slo for B2u- (o) 
and A2u-symmetric (x) modes of the same frequency. Because the cuq is far from the 
characteristic frequency of the sign-change (Fig. 9(b)), the onset is very abrupt. 

Fig. 12: (a) Phonon linewidth broadening as a function of T for hypothetical Big- 
(o) and Aig-symmetric (x) modes of frequency ujq = 0.2 J, at g = for the YBCO-like 
system. Because ujo is close to the frequency of the peak (Fig. 10(a)), a significant 
broadening is observed, (b) ST for B2u- (o) and A2„-symmetric (x) modes the same 
frequency. Because uq is so far from the characteristic frequency of the sign-change 
(Fig. 10(b)), the imaginary part is negligible. 
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